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We consider the angular momentum of a harmonically trapped, noninteracting Fermi gas subject 
to either rotation or to an artificial gauge field. The angular momentum of the gas is shown to 
display oscillations as a function of the particle number or chemical potential. This phenomenon is 
analogous to the de Haas - van Alphen oscillations of the magnetization in the solid-state context. 
However, key differences exist between the solid-state and ultracold atomic gases that we point 
out and analyze. We explore the dependence of the visibility of these oscillations on the physical 
parameters and propose two experimental protocols for their observation. Due to the very strong 
dependence of the amplitude of the oscillations on temperature, we propose their use as a sensitive 
thermometer for Fermi gases in the low temperature regime. 

PACS numbers: 05.30.-d, 03.75.Ss,05.30.Fk,67.85.-d 

Keywords: Quantum oscillations, ultracold Fermi gases, rotating atomic gases, artificial gauge fields 



I. INTRODUCTION 

When subject to an applied magnetic field, the elec- 
tron gas of a metal displays an oscillatory dependence on 
physical observables as a function of field strength. Gen- 
erally known as 'quantum oscillations', these effects can 
be observed for example for the magnetization (de Haas- 
van Alphen oscillations), the resistivity (Shubnikov-de 
Haas oscillations), the Hall resistance, or the specific 
heat. Quantum oscillations are a macroscopic manifesta- 
tion of the Landau level quantization of the spectrum of 
electrons in the presence of a magnetic field. They were 
predicted by Landau [1] and observed by de Haas and 
van Alphen [2] in 1930. Since then, quantum oscillations 
have become an invaluable tool in characterizing the 
electronic states of metals and semi- metals with large 
mean-free path. The measurement of the oscillations 
allows in particular for a quantitative determination of 
Fermi surfaces and of the effective masses of quasiparti- 
cles. A recent illustration is the remarkable observation 
of quantum oscillations in underdoped cuprates in large 
magnetic fields [3]. This measurement has triggered 
many discussions about the reconstruction of the Fermi 
surface in these systems and the possible relevance of 
Fermi liquid descriptions. In this article, we explore 
the possibility of observing similar quantum oscillation 
effects for ultracold gases of fermionic atoms, either by 
putting the gas in rotation or by using artificial gauge 
fields. We focus mostly on the total angular momentum 
of the g key observable, and identify the regimes 

in which oscillations of this quantity as a function of e.g. 
atom number are pronounced enough to be observable 
experimentally. We emphasize the similarities, but also 
the key differences, with the solid-state context. 
Rotation has proven to be an efficient way to mimic 
the effect of a magnetic field for ultracold gases of 
neutral atoms. Over the last years, strong experimental 



and theoretical effort [4-6] has been devoted to the 
understanding of rotating cold atomic gases. Experi- 
mental achievements include notably the observation 
of vortex nucleation in Bose-Einstein condensates [7] 
and Fermi gases [8], which has been used as a proof 
of their superfluidity. On the theory side, rotating 
gases are expected to be highly controllable test ta- 
bles for strongly correlated phases of matter under 
the influence of a magnetic fleld [4-6], such as the 
fractional quantum Hall regime. Recently, the experi- 
mental realization of artificial gauge fields with Raman 
beams [9, 10[ opened another promising way towards 
the generation of artificial magnetic fields for cold atoms. 

In this paper, we analyze the infiuence of rotation or 
of an artificial magnetic field on a noninteracting Fermi 
gas, in a cylindrically symmetric harmonic trap. We fo- 
cus mostly on the angular momentum of the gas with 
respect to the rotation axis. The angular momentum 
of a quantum gas differs from its classical counterpart. 
While the former is defined from the change in energy 
when varying rotation frequency, the latter is given by 
the moment of inertia of the gas considered as a rigid 
body and hence is related to its mean-square radius in 
the plane perpendicular to the rotation axis. We compare 
the quantum angular momentum to its classical counter- 
part by introducing their ratio. This quantity has already 
been studied in the context of rotating Bose-Einstein con- 
densates. When Bose-Einstein condensation occurs, this 
ratio deviates from unity [11[. It has also been shown 
to be proportional to the difference of frequencies of the 
m — 2 quadrupolar oscillation modes [7, 12, 13[. In the 
case of superfluid helium, this ratio can give access to the 
ratio of the normal and superfluid fractions [14, 15[. 

In the fermionic gas considered in this article, the ratio 
of the quantum angular momentum to its classical coun- 
terpart tends to unity for very large atom number. We 
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show here that for intermediate values of the atom num- 
ber and at low enough temperature, this ratio displays 
pronounced quantum oscillations as a function of atom 
number or chemical potential (Sec. II, Sec. III). Consid- 
ering both the slow and fast rotation regime, we provide 
(Sec. IV) a simple qualitative picture of these oscillations 
and emphasize the similarities but also key differences 
with the de Haas-van Alphen oscillations in the solid- 
state context [16]. We perform a detailed analytical and 
numerical analysis of the dependence of the amplitude 
and period of these oscillations on temperature, rotation 
frequency, atom number or chemical potential and as- 
pect ratio of the trap (Sec. V). We emphasize in partic- 
ular that this oscillatory behavior, due to its sensitivity 
to temperature, may be used as an accurate thermome- 
ter in the low (T/Tp < 0.1) temperature regime. Finally, 
we propose (Sec. VI) two different experimental protocols 
that aim at measuring the angular momentum of the gas, 
and thus to reveal its oscillating behavior with respect to 
particle number. 



II. MODEL AND OBSERVABLES 

A. General framework 

We consider a noninteracting Fermi gas in a cylindri- 
cally symmetric harmonic trap (with the radial and axial 
trapping frequencies being loq and cj^, respectively), ro- 
tating at an angular frequency uj around the symmetry 
axis Oz. The Hamiltonian in the rotating frame has the 
following form [4-6]: 



n 



2M 



1 



2 „2 



(1) 



where p is the momentum, AI is the atomic mass and 
is the angular momentum with respect to the rotation 
axis. The Hamiltonian (1) is an analogue of the Hamil- 
tonian describing an electron moving in a magnetic field 
in the solid state context. The angular momentum is the 
analogue of the magnetization, and the angular frequency 
plays the role of the magnetic field. However, note that 
typically no external trapping potential is present in the 
solid state context. We will see that this can lead to quite 
different behavior. 

The Hamiltonian (1) gives rise to the following spec- 
trum [4, 17]: 

= fiuj+ (n+ + l/2)+hLj^ (n_ + l/2)+huj2 [n^ + 1/2) , 

(2) 

where LiJ± = ujq^w and v — (7i+, n_ , n^) with the nonneg- 
ative integers n+,n_,n^. Interestingly, (2) is identical to 
the spectrum of a noninteracting Fermi gas in a fully 
anisotropic harmonic trap, with trapping frequencies ljj± 
and Uz- 

The gas will be described at fixed inverse temperature 
/? = -jT^, chemical potential /i, and angular frequency lo 



in the grand canonical ensemble through the thermody- 
namic potential f2 : 

n = n{p, fi, a,az)^ -ksT log [l + e'^^""^'')] , (3) 

M 

We defined the following dimensionless quantities : 

a = — , and the aspect ratio of the trap az — — , 
ujQ Wo 

(4) 

which together with will be the relevant parameters 
of the problem. 



B. Angular momentum 

The angular momentum {Lz) of the gas is given by 
the derivative of the thermodynamic potential (3) with 
respect to the angular frequency cj : 



{Lz) 



du! 



(5) 



Here, it reduces to: 



^ h{{n^) - {n+)) . 



(6) 



In the first equality, /(e) = jq:^ 
function at temperature T. 



denotes the Fermi 



The angular momentum of a classical gas is propor- 
tional, in a rigid body-like motion, to the classical mo- 
ment of inertia, which is proportional to its square ex- 
tension around the rotation axis z (cf. appendix C) : 



(7) 



It should be mentioned that this expression remains true 
even for interacting particles, as shown in appendix C. 

In order to derive a quantum analogue of the classi- 
cal expression, we use that the square extension is given 
by the derivative of 51 (cf. eq. (3)) with respect to the 
squared radial trapping frequency • 

(8) 

This leads to the following expression for the quantum 
analogue of the classical angular momentum : 

(L,) = /ia ^ /(e, - ;i) {n+ + n_ + 1) (9) 

= ha + (n^) + N) , (10) 

where N = li^v^lA is the total number of particles 

of the gas. 
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In order to quantify the deviation of (6) from its clas- 
sical counterpart (10), let us define the ratio 



R 



(L,)M) a{{n^) + {n+)+N)' 



(11) 



which depends on temperature, chemical potential (or 
equivalently on particle number), rotation, and axial con- 
finement. In the limit of high temperature and large 
sample size, the gas behaves classically and the angular 
momentum is well described by its classical expression 
(10), and consequently R= \. Deviations from this clas- 
sical analogue will be interpreted as manifestations of the 
quantum behavior of the gas. The analysis of these de- 
viations is the main subject of this work. 



III. OSCILLATIONS IN THE ANGULAR 
MOMENTUM 

Fig. 1 depicts the evolution of R with the atom number 
for a rotating fermion gas. 

In this plot, one sees that the difference between the 
angular momentum of the gas (5) and its classical ana- 
logue (7) versus the particle number N is twofold : 

• A monotonous deviation Rmom given in eq. (17) 
[18] 

• An oscillating contribution Rose with N which is 
very pronounced at low particle number . 

In order to understand these two different contri- 
butions, we step back and investigate the thermody- 
namic potential $7 (3) from which the physical quantities 
will be extracted. Appendix A shows how the expres- 
sion (3), via Poisson summation, splits into two parts 
n = + The first one denotes the non os- 

cillating contribution with the chemical potential and is 
given by : 



-Li4[— e 



/3(p-£o)l 



(12) 



The second part is responsible for the oscillations with 
the chemical potential and is determined by 



2kBT'^{ > with 

[0,+ooP 



(13) 



(14) 



In the previous equations, Li„ is the n-th poly logarithm 
function [19], k = {k+, k^) is the vector conjugated 
to u = (u+, u-,Uz) with fc-f, fc_ and kz positive integers, 
eo = ^0 + \ f^z is the ground state energy and = 

The splitting of the thermodynamic potential into an 
oscillating and a monotonous contribution with /i carries 
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Figure 1. The ratio R, for = 0.05, a = 0.01, and T /Tf = 
0.03 (in blue). The black dashed line is the monotonous con- 
tribution, and the red dotted line is the oscillating part. At 
high particle number, the oscillating contribution is negligi- 
ble, and R coincides with its monotonous part (17), as they 
both go to unity. 



over to many observables. Typically, these oscillations 
stem from 'shell filling' effects as we will show in the 
next section. They have been predicted for different 
observables as e.g. the radius of the cloud [20-24]. 
However, the relative strength of both contributions 
can vary very much with the observable considered. 
We will show that the oscillating contribution is much 
more pronounced in the angular momentum {L^) of the 
gas with respect to the rotation axis then in previously 
considered observables [20-24]. 

From the splitting of the thermodynamic potential into 
the two contributions, the analogue separation of both 
{Lz) and its classical analogue follows : 



R = 



(L,)(«°) + (L,)(°-) 



— Rri 



Rn 



The monotonous contribution is given by : 



Rn 



(L^)(ci,no) 



(15) 



(16) 



2Li4 



j^(M-eo)j 



2Li4[-e/3(^-^o)] + (/3fia;o)4(l - a^^azN ' 

(17) 

The expression (17) goes to unity for large samples as 
depicted in Fig. 1 which means that the classical value 
of the angular momentum is reached. The variations of 
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the ratio R with respect to N, a, and az stem both 
from a direct dependence on these parameters as ev- 
ident in (17), but also from the indirect dependence 
via the chemical potential. For example, at zero tem- 
perature the chemical potential can be expressed by 

^(T = OK) = hujo {6{l ~ a^)azNy^^. Using this ex- 
pression and developing the Li„ functions in (17), one 
sees that the strongest dependence is seen at low tem- 
peratures, where the expression simpUfies to: 



1 



l + 2(l-a2)2/3(6a^iV)-i/3 



(18) 



The oscillating part of R (the red dotted line in Fig. 1) 
is recovered by expanding the denominator in (16) in 
power series. It is of the general form : 



Ros 



23? 



^i?(fc)e'*"^ 




(19) 



where R{k) contains contributions from (L^ 

The expression (19) states clearly that this contribu- 
tion has a periodic dependence on the chemical potential, 
and therefore on the particle number, and that the char- 
acteristic frequencies, that are the oscillation periods, are 
Ltj± and ojz- The next section will shed some light on 
this oscillating behavior, on the basis of qualitative argu- 
ments. 



IV. QUALITATIVE INTERPRETATION 

This section provides an intuitive interpretation of the 
oscillatory behavior of the angular momentum at low 
temperature. We shall first review a simple physical pic- 
ture of quantum oscillations in a two-dimensional elec- 
tron gas, and then point out the differences and similar- 
ities in the case of rotating cold atomic gases. We derive 
in particular simple analytical expressions for the angu- 
lar momentum in the case of a two-dimensional system 
— >■ oo) in the regime of fast rotation w < wq- 



A. Quantum oscillations in a two-dimensional 
electron gas: a simple picture 

Let us first consider the case of electrons in a strictly 
two-dimensional electron gas (2DEG). Quantum oscil- 
lations have been reported in such systems, see e.g. 
Ref. [25]. In this case, such effects can be understood 
from a simple calculation which we now recall [16]. 

In the presence of a magnetic field B, the energy lev- 
els of the 2DEG are Landau levels with single-particle 
energies En = (n + l/2)hLUc, with rt > an integer and 
(jjc = ^ the cyclotron frequency (m is the mass of the 
electron and — e its charge). Each Landau level is macro- 
scopically degenerate, with the degeneracy given by 



where $ = BS is the magnetic flux over the whole sam- 
ple, with B the amplitude of the magnetic field and S the 
surface of the sample, and $o = ^ is the flux quantum. 

Let us consider the system at zero temperature with N 
electrons, in the case where the first p Landau levels are 
completely filled (each one with ^/^q electrons) and the 
(p-l-l)-th one is partially filled with N —p^/^o electrons. 
The energy of the system E = T,lJo{n + 1/2) hujc^/<^Q + 
{N - p^/^o)hujc{p + 1/2) reads: 



huJc 2 

bj ~ p 

2 $0 



N -p^]{p+l/2)tiuj,. (20) 
^0 , 



Taking a derivative, one obtains the magnetization of the 
electron gas for < N < {p + 1)^: 



M = -g^=M2p+l) 



$0 



A^BP^ (21) 



with fJ-B — ^ the Bohr magneton. This expression 
shows that the magnetization jumps discontinuously by 
a macroscopic amount 2N each time a new Landau 
level starts to be filled: when exactly p Landau levels are 
filled the magnetization per electron is M/Nfig = — 1 
(Landau diamagnetism) , and it jumps to AI/Nfis — +1 
as the {p+l)-th level starts to be filled, reaching smoothly 
M/NjiB = — 1 again when this new level is completely 
filled. This behavior is depicted on Fig. 2, which dis- 
plays M/NjiB as a function of N^q/^. These abrupt 
jumps are due to the macroscopic degeneracy of the Lan- 
dau levels. The magnetization per electron at T = is 
an approximately periodic function of N^q/^, that is of 
particle number and inverse magnetic field. 

While this picture applies to 2DEGs, the electron dis- 
persion along the z-axis is of course not negligible in 
most solid-state materials and must be taken into account 
(leading e.g. to the semi-classical picture of quantum os- 
cillations in terms of extremal closed orbits on the Fermi 
surface) . However, the simple picture above captures the 
essence of the de Haas - van Alphen effect. 



B. Physical picture for a trapped rotating Fermi 

gas 

1. Geometrical representation 

We now investigate the case of a rotating trapped 
atomic gas of fermions at T = 0, focusing again for sim- 
plicity on a purely two-dimensional gas (strong trans- 
verse confinement ujz — >■ oo). The main difference with 
the case of a 2DEG is that the single-particle spectrum 
is now that of a two-dimensional anisotropic harmonic 
oscillator, with two different frequencies ij± = ujq ± cj. 
Hence, energy levels are non-degenerate, in contrast to 
the macroscopic degeneracy of the Landau levels typi- 
cally present in the solids. 
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the particles with respect to the order in which states are 
filled. 

Increasing the chemical potential /i will lead to the fill- 
ing of new states. In the situation sketched in Fig. 3 the 
difi'erence between the two different chemical potentials 
corresponds to exactly adding one particle, namely par- 
ticle number 16. Since this particle has negative angular 
momentum, it will lead to a decrease of the total angular 
momentum (L^) (cf. Fig. 3). 

This simple picture shows very clearly the origin of the 
oscillations in the angular momentum. Roughly speak- 
ing, depending if during a shift of the chemical potential 
more particles are added in the upper diagonal (resp. 
lower) diagonal part, the angular momentum will in- 
crease (resp. decrease). For particles on the diagonal 
the angular momentum does not change. 

2. Simple picture at fast rotation 



Figure 2. Dashed curve (red Hne, left scale): magnetiza- 
tion per particle M/NfiB of a two-dimensional electron gas 
at zero-temperature, as a function of N^o/^, with TV the 
number of electrons and $/$o the magnetic flux in units of 
the flux quantum = h/e. A discontinuous jump is ob- 
served each time a Landau level is completely flUed. Plain 
curve (blue line, right scale): dependence of the angular 
momentum L^/N'^h on Nlo-/lo+ for a zero-temperature ro- 
tating trapped fermionic gas, in the limit of fast rotation. 
The angular momentum displays cusps at the discrete val- 
ues Nui-/ujj^ — p(p + l)/2 = 1, 3, 6, 10, • ■ • , corresponding to 
the complete flUing of an integer number of n+ levels. At 
these points, the derivative of the angular momentum ver- 
sus Nuj-/uj+ displays discontinuous jumps. This illustrates 
the similarities, but also the differences between the solid- 
state and ultra-cold fermionic gas contexts, regarding the phe- 
nomenon of de Haas - van Alphen quantum oscillations. 



Fig. 3 helps to understand how this reflects on the 
angular momentum. This picture displays the single- 
particle states on a two-dimensional grid labeled by 
(n-)_ , n_ ) . Because the energy levels (counted from the 
zero-point energy Hujo) read huj-\-n^ + hu^n^, the value 
of the energy of a given state is obtained geometrically by 
projecting the representative point (n+,n_) onto an axis 
tilted by an angle tan0 — uj-/uj^. Since each state can 
be filled by at most one fermion, the ground-state con- 
figuration for a given chemical potential (Fermi level) /.t 
is obtained by occupying all states whose projection onto 
this axis is lower than /i. The angular momentum of a 
single-particle state, on the other hand, being given by 
iz/h = fi- — n+ is obtained by projecting the point onto 
the anti diagonal of slope —1 (Fig. 3). This means that 
particles exactly on the diagonal will not contribute to 
the angular momentum, whereas particles close to the 
vertical (ti-) and horizontal («.+ ) axis will contribute 
most, with a positive and negative sign, respectively. The 
evolution of is shown in the inset of Fig. 3, ordering 



For the case of fast rotation (w « wq, so that a;+ >> 
a;_) this reasoning can be made more quantitative, and 
a simple calculation in the spirit of the one on 2DEGs 
above can be done. In this case, an analogue of the Lan- 
dau levels are the successive n+ levels. Further, the factor 
w+ /a;_ relates the number of n_ levels with respect to the 
number of n_|_ levels analogous to $/5'o, so that fast rota- 
tion corresponds to large magnetic fields. However, a cru- 
cial difference is that each n+ level is not macroscopically 
degenerate, since aj_ is small but non-zero. The number 
of available n_ states is nevertheless much larger than 
the number of available states. Therefore, we can as- 
sume a continuum of levels, in the n_ direction, while 
treating the n+ levels are discrete. Assume that at a 
certain chemical potential (particle number), the largest 
allowed n+ quantum number is given by p. Then the 
largest allowed 7i_ which corresponds to a certain fixed 
n+ is given by a;_(n+) — uj^/uj^(jp + 5^{N) — n_|_), where 
5+{N) = {N ~lo+/uj^p{p+1)/2)/{lu+/lu^{p+1)) is the 
distance of crossing of the chemical potential line and the 
axis from p, as depicted on Fig. 3. 

In this picture the angular momentum can be calcu- 
lated as: 



E 



n+=0 



:-("+) 



(n_ — n^) 



(22) 
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p{l+p) 
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In the limit of very fast rotation a;-|_ ^ this can 
be cast in the simpler form, valid for p{p + l)/2 < 
iVa;_/w+ < {p+l){p + 2)/2: 
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Figure 3. Sketch explaining the origin of the quantum os- 
cillations for a two-dimensional gas at zero temperature (the 
case of low confinement is similar. In this case, each point 
of the graph will account for the particles that populate the 
tJz levels allowed by the energy constraint). Each level is la- 
beled by its integer coordinate (n+, n_) and a projection onto 
the -& axis with slope uj+/uj- gives the single particle energy 
£ — eo. The occupation of the levels is fixed by the chemi- 
cal potential (green line). Green circles represent filled levels 
at zero temperature (satisfying the energy constraint), and 
the green situation depicted corresponds to the ground state 
of the system with N=15 particles. The angular momentum 
corresponds to the projection onto the axis labeled by l^/h as 
indicated by the dashed-dotted line. The two different colors 
for the axis and the corresponding shaded area indicate differ- 
ent sign of i^. The particles have been attributed a number 
(inside the green circle) corresponding to their order of ap- 
pearance as the chemical potential is increased. Inset: The 
inset shows the evolution of the total angular momentum of 
the gas (in blue) with the number of particles, and the an- 
gular momentum of each particle (orange). The arrows point 
the total angular momentum for A'^ = 16, and the contribu- 
tion from the supplementary particle added to the system by 
changing the chemical potential from the green to the red line. 



generacy is lifted due to the presence of the w_ oscillator. 



3. Oscillations at slow rotation 

From this picture, and from the geometrical represen- 
tation above, we expect that stronger oscillations are ex- 
pected to occur in situations in which the number of occu- 
pied + and — are comparable, i.e. for slow rotation. Fur- 
thermore, leaving the assumption of a two-dimensional 
system, we see that many particles can be added in the 
uiz levels between successive uj± levels if uj^ is small, thus 
leading to oscillations occurring for larger systems. 

This simple scheme also explains why the oscillations 
are much less pronounced for the classical angular mo- 
mentum (rigid-body moment of inertia). The u!± levels 
contribute to this quantity with the same sign (Eq. (10)). 
Hence, only shell-filling effects for very small systems 
are expected, as previously discussed for the radius of 
a mesoscopic gas confined to an anisotropic trap [20- 
24]. These oscillations, e.g. in the extension of the cloud 
are very hard to resolve experimentally and are typically 
only visible for very small atom numbers. Fig. 4 com- 
pares, for weak confinement, the effects in and in 
its classical analogue (extension of the gas). It is seen 
that in contrast to the latter, Lz displays pronounced 
oscillations as a function of the number of particles for 
rather large systems. This is because the angular mo- 
mentum (Lz) involves the difference of n+ and n_. This 
also shows clearly that the main oscillatory contribution 
in (11) stems from {LzY°"'\ 

In this section, we have gained a helpful intuitive pic- 
ture of the occurrence of oscillations in the angular mo- 
mentum. The following section explores in more details 
the influence of the physical parameters of the system on 
the deviations of the angular momentum from its classi- 
cal counterpart, both analytically and numerically. 



INFLUENCE OF PHYSICAL PARAMETERS 
ON QUANTUM OSCILLATIONS 



which is to be compared to the form (21) for a 2DEG. In 
contrast to this case, it is easily checked that this formula 
is actually a continuous function of Nlo^/ujj^- (Fig. 2). 
At the special filling factors Nuj_/uj+ — p{p + l)/2 
corresponding to starting to fill the p-th level of the 
a;+ oscillator, Lz does not display a discontinuous 
jump but rather a cusp, at which only the derivative 
of Lz with respect to Nuj_/uj^ has a discontinuity. 
This is clearly visible on Fig. 2. This difference in the 
behavior of (Lz) and the magnetization M of an electron 
gas, which are analogous quantities, stems from the 
macroscopic degeneracy of the Landau level in the case 
of electrons. In the case of the rotating gas, the de- 



In the forthcoming section, we will essentially focus on 
the oscillating contribution to the ratio of the quantum 
and classical angular momentum R. We discuss the im- 
pact of temperature, axial confinement and angular fre- 
quency on the amplitude and period of the oscillations, 
through the dimensionless quantities Uz and a, re- 
spectively. In an experiment, typically the atom number 
can be measured more easily than the chemical potential. 
Therefore, we monitor the evolution of physical quanti- 
ties as functions of the atom number rather than the 
chemical potential. Since the amplitude of the oscilla- 
tions will vary considerably with the system parameters 
we will point out the regimes in which the oscillations are 
best visible at reasonably large particle numbers. 
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Figure 4. Angular momentum and its classical counterpart as 
a function of iV, for = 0.05, T/Tf = 0.01 and a = 0.05. 
The shell filling effect in the angular momentum is magnified 
due to the difference between n+ and n_ in its expression, 
and results in an oscillating behavior. The inset shows the 
same two quantities, at low particle number. The shell filling 
effects are visible in the classical angular momentum only at 
low particle number. 



Figure 5. R vs UzN , at low temperature T/Tf = 0.03 and 
a = 0.05, for = 0.01, = 0.02, = 2 and = 10. Even 
though at low particle numbers, the amplitude of the oscil- 
lation is strongest for the intermediate aspect ratio ct^ = 2, 
the oscillations remain visible for larger particle numbers in 
the weak confinement case. This indicates that a compro- 
mise between these two feature might be a favorable choice of 
parameter to observe the oscillation of angular momentum. 



A. The oscillation amplitude 

In this section, we focus on the influence of the system's 
parameters on the amplitude of the quantum oscillations. 



1. Axial confinement and temperature 

Fig. 5 compares the angular momentum ratio i? as a 
function of particle number for different values of the as- 
pect ratio Uz of the trapping, and for identical values of 
angular frequency and temperature. For all parameters, 
the monotonous contribution leads to a rise of the ratio 
R with increasing particle number towards the value one 
signaling a classical behavior of the angular momentum. 
On top of this rise, clear oscillations are evident at low 
particle number. The strong damping of the oscillations 
at large particle number is expected from the physical 
intuition given by Fig. 3. The occupation of a new level 
becomes less important due to the large number of pos- 
sible configurations. Also the denominator in (11) has 
larger values, and consequently a damping associated to 
the increasing particle number is visible. 

The confinement along the ^-direction controls how 
many levels are already at play before new -t-, — levels 
contribute (cf. Fig. 3). Therefore, the sharpness of the 



oscillations is strongly influenced by the aspect ratio of 
the trap as evident in Fig. 5. The largest amplitude of the 
oscillations is found at an almost isotropic confinement 
and low particle number N . However, at lower aspect 
ratio the oscillations survive till larger particle numbers, 
which might be a favorable situation to observe quantum 
oscillations in experiments. 

In the low case, the oscillations reach a universal 
regime, where the curves coincide as functions of UzN , 
as depicted on Fig. 5. This dependence on a^N can be 
understood transforming the sum over in (3), into an 
integral, which gives for the thermodynamic potential : 

^^{kBTY_ r_g/3(M-eo-ft^+n+-R^-«-)l ^ (25) 

tu^Jz ^ 

The only dependences on cj^ are in the prefactor and in 
Eg. However, in the weak confinement regime, ui^ ^ wqj 
and consequently the zero-point energy is dominated by 
Eq ~ hujQ. Taking the ratio R, the prefactor drops out 
in the expression (11), and one is left with a function 
of /i oc {a^NY^^ only. This universal regime strongly 
resembles the van Alphen oscillations in a solid. Due to 
its large extensions in a solid typically the electrons are 
assumed to move freely in the direction of the magnetic 
field. 

The previous reasoning is valid at low temperature. At 
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Figure 6. (a)- The ratio R vs particle number, at low rotation frequency a = 0.05 for T/Tf = 0.03, T/Tf = 0.05, T/Tf = 0.07 
and T/Tf = 0.1. Panel (a) shows a weak aspect ratio of the confinement = 0.02 and panel (b) an intermediate aspect 
ratio az ~ 2. In both cases, the low temperature curve displays the most pronounced oscillations, and the oscillating behavior 
remains visible for larger samples. The stronger confinement leads to larger oscillation amplitudes, but the oscillations disappear 
for lower particle numbers. The insets shows the Fourier transform of R wrt fi for the corresponding parameters. It shows that 
the amplitude of the oscillation modes decreases with the temperature, but the oscillation period remains unchanged. 




Figure 7. Amplitude of the oscillation of the ratio R vs T/Tf, 
for a = 0.05 and for a^N = 11.5 (blue) and q^A'^ = 4.5 
(black) with respect to the monotonous contribution (16). As 
expected, at fixed confinement, the oscillation becomes less 
pronounced for bigger particle numbers. The strong damping 
effect due to the temperature is clearly visible in both cases. 



larger temperature, the sharpness of the Fermi surface is 
lost which will lead to less pronounced oscillations. Fig- 
ures 6- (a) and 6-(b) compare the oscillations for different 



temperatures and fixed rotation, in the weak and strong 
confinement cases, respectively. A strong damping with 
increasing temperature of the amplitude of the oscilla- 
tions is evident. The strong dependence of the oscilla- 
tions on the sharpness of the Fermi surface plays an im- 
portant role in the observability of quantum oscillations, 
as known from the solid-state context [16], and pointed 
out for cold atoms in [20, 22-24]. 

Fig. 7 presents the amplitude of the oscillations as 
a function of temperature, at fixed rotation and axial 
confinement. A strong increase at low temperatures 
which seems almost linear occurs. This strong tem- 
perature dependence of the amplitude of the oscillat- 
ing contribution can be understood from a Poisson sum- 
mation analysis as presented in appendix A, applied to 
(Lz) and (Lz)^'^''. By doing so, the temperature de- 
pendence only occurs in the sum in form of the expres- 



( sinh 



TT Kb 



This term 



causes the strong temperature damping of the oscillating 
behavior of thermodynamic quantities. Let us note, that 
it also shows that the temperature has no influence on 
the oscillation period as we discuss below. 

The strong dependence of the oscillation amplitude on 
temperature offers the possibility to use the oscillations 
as a thermometer in the low temperature regime. This 
low temperature regime (T/Tp < 0.07) is experimentally 
challenging, but has nevertheless been realized in a dif- 
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ferent context, with rotating [8] or non rotating [26, 27] 
apparatus. A thermometer in this regime would be very 
important, since commonly used thermometry as the mo- 
mentum distribution after time-of-flight imaging becomes 
less sensitive to temperature changes in this regime. 



2. Rotation frequency 

In the solid-state case, quantum oscillations are most 
easily seen with high magnetic fields, which translate here 
into fast angular frequency. This statement finds its ori 
gin in the criterion that the cyclotron frequency a;, 
must be bigger than ksT. 
Fig. 8 depicts the evolution of R with the particle num- 
ber for different values of a at constant T jTp and axial 
confinement. In contrast to the solid state case the os- 
cillations here are more pronounced when the rotation 
of the gas is slow, as already discussed qualitatively in 
section IV. This can be related to the occurrence of the 



_B 
M 



two characteristic frequencies 



Indeed, a should 



not be too high such that remains bigger than fc^T. 
This is expected from the Poisson summation analysis : 
lower values of hw^ (corresponding to higher values of 
a) tend to increase the temperature damping factor (see 
Appendix A) and thereby decrease strongly the ampli- 
tude of the oscillations. Thus a strong dependence on 
the frequency w is found. 

Although a must be low compared to 1, the rotation 
must not be too slow to be observable. 



B. (Pseudo)Period of oscillation 

In this subsection we focus on the effect of experimen- 
tal parameters on the pseudo period of the oscillations. 
From the intuitive picture described in section IV, one 
expects that the period of the oscillations does mainly 
depend on the energy level spacing given by huj± 2- 

The strong dependence on the two fundamental 'fre- 



quencies 



and 



t, f-!— is confirmed by our numerical 

investigations reported in Fig. 9. This figure displays the 
modulus of the Fourier transform of the ratio R (11) for 
several angular frequencies at zero temperature (upper 
panel) and for T/Tp = 0.03 (lower panel). The peaks 
that develop around integer values of xhujo (the variable 
conjugated to n) show that two fundamental 'frequencies' 
and appear (pointed out by arrows in Fig. 9). 
Increasing the frequency w the peaks become more and 
more separated. Additionally to these main peaks, a peak 
at shows up here well separated from the other peaks 
(see inset in the upper panel of Fig. 9) . The interplay of 
these three different frequencies is also well seen in the 
quantum oscillation shown versus the particle number in 
Fig. 8. The influence of u}± can in particular be identi- 
fied in Fig. 8 comparing a = 0.1 and 0.2. A clear shift of 
the main oscillation period occurs. Additionally to the 




500 



1000 




a=0.15 

a=0.2 



500 



1000 



Figure 8. R vs. particle number, at zero temperature (up- 
per panel) and T/Tp = 0.03 (lower panel) and weak trapping 
ftz = 0.075, for several angular frequencies. Slow rotation 
favors the appearance of quantum oscillations. At zero tem- 
perature, fast oscillations due to the contribution of the uiz 
levels appear. 



frequencies uj± , the trapping frequency causes smaller 
oscillations (Fig. 8). Chosing a larger value of the aspect 
ratio of the trap, will infiuence considerably the oscilla- 
tion period as shown in Fig. 5. 

Most of the detailed substructures of the oscillation 
discussed above are washed out at higher temperature, 
since the temperature causes the sharp peaks in the 
Fourier transform to broaden and only the main fre- 
quency determined by ujq survives, as shown in Fig. 8 
and 9. 

One major difference compared to the solid state con- 
text is the existence of the three characteristic frequencies 
h{ujl±uj) Tuir ™ ^^'^ oscillations, which is an effect 
of the radial and axial trapping potential, respectively. 
This occurrence is responsible for an intrinsic limitation 
on angular frequency, which is forced to be lower than 
Wq so that a;„ remains positive. This limitation is also 
at the origin of the difficulties encountered to generate 
artificial magnetic fields strong enough to observe the 
quantum Hall effects in rotating gases [4-6]. One possi- 
bility to circumvent this problem is to add a quartic term 
to the trapping potential, which has been implemented 
for example in [28]. This intrinsic limitation is also a mo- 
tivation for the implementation of artificial gauge fields, 
which form the basis of one of the protocols described in 
section VI. 



10 



0.03 - 




0.02 



r 




1 

2 


1 

3 


1" 
4 






1 


1 






xhuo 



Figure 9. Fourier transform of R vs xfaoo at zero tempera- 
ture (upper panel) and T/Tf = 0.03 (lower panel) , for fixed 
Qfz = 0.075. The successive harmonics of the fundamental fre- 
quencies are visible. The arrows indicate the position of the 
two characteristics 'frequencies' jr^. Inset : Fourier trans- 
form of 71 at zero temperature showing the peak correspond- 
ing to -r—, indicated by an arrow. 



VI. MEASUREMENT OF THE ANGULAR 
MOMENTUM 

In section V, we have shown that the quantum oscil- 
lations of the angular momentum are best visible at low 
temperatures {T/Tp < 0.07) and slow rotation (a < 0.2). 
A compromise between the particle number and the axial 
confinement has to be found. Even if hard to fulfill, these 
requirements are nevertheless on the reach of current ex- 
perimental capacities. The height of the amplitude of the 
oscillations will be a very sensitive probe to the temper- 
atures of the gas. 

In order to measure R, one has to measure two observ- 
ables : the square radius of the gas, and its angular mo- 
mentum. Measurements of the density of the cloud and 
thus of its square radius can be performed using spatially 
resolved in-situ imaging techniques which are nowadays 
standard. For the observation of (i^) two different ex- 
perimental protocols (Fig. 10) are useful. Since these are 
not commonly used, we briefiy summarize them in the 
following section : 

• Exciting the quadropole modes of the rotating gas 
and measuring the induced precession angle [7, 29] , 
depicted in Fig. lO-(A) 

• Measuring the occupation of internal states of the 
atoms which are connected to the angular momen- 



tum for light induced vector fields [30, 31], shown 
in Fig. lO-(B) 

The first method has been successfully implemented 
with rotating Bose-Einstein condensates ]7, 29]. The un- 
derlying idea, as depicted in Fig. lO-(A), is to imprint an 
anisotropic perturbation of the form: 



y 



(26) 



for a time 5t shorter than the rotation period lo^^ ]32]. 
This perturbation will cause an excitation of quadrupo- 
lar oscillation modes, which will lead to a precession of 
the cloud. The precession angle Lp of the main axis of 
the gas (see Fig. 10- (A)) is related to the geometrical 
properties of the gas : 



tan2(/3 



2{xy) 



(.t2) - 



(27) 



After the short excitation these properties change with 
time and can be computed perturbatively in the param- 
eter e 



Vo&t. 

h 



{x^) - {y^){t) ^ ^-^rlt 



eh 
M2 
2he 



(28) 
(29) 



where r§ is the radial size of the gas before the pertur- 
bation is applied. Assuming that is small leads to the 
relation: 



LP 



2Mrl 



t. 



(30) 



Thus, following the evolution of this precession angle with 
time provides a direct measurement of 

The second technique to measure the angular mo- 
mentum is described in Ref. [30, 31], making use of 
artificial gauge fields to generate an artificial slow 
rotation. More precisely, one has to shine on the gas 
two CO propagating Laguerre-Gauss beams carrying 
different angular momenta ^ and i -\- A^, as indicated in 
Fig. lO-(B). By two-photon Raman transition, the atoms 
become dressed by these beams and the dispersion of 
the lowest band is displaced by an effective shift to a 
non-zero angular momentum that is proportional to A£. 
Since the generation of the vector field by this gauge 
field induces a coupling between the external (angular 
momentum) and internal (hyperfine) states of the atoms, 
the population of the internal states can be used to infer 
the angular momentum of the cloud. 

Both these methods can be used as a thermometer of 
the fermionic cloud at low temperature, as suggested in 
section V. By comparing directly the measured oscilla- 
tion pattern with the exact expression derived from the 
thermodynamic potential (3), knowing all the other pa- 
rameters allows to determine the temperature of the gas. 
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Figure 10. Sketch of the two different protocols for measur- 
ing the angular momentum of the fermion gas (see text for 
description): (a) Measurement of the precession angle. The 
stirring beam (in green) is also used to create the anisotropic 
perturbation, which provokes the appearance of quadrupole 
modes, and are responsible for a precession angle ip. Observ- 
ing this angle provides a direct measurement of the angular 
momentum of the gas . (b) Generation of artificial rotation 
with lasers with spectroscopic measurement of the angular 
momentum. The two co propagating beams carry different 
angular momenta, and the effective rotation is proportional 
to the angular momentum difference. The angular momentum 
is then measured by spectroscopy as suggested in [30, 31]. 



Furthermore, as Fig. 7 has shown, the amplitude of the 
oscillatory phenomenon is expected to depend strongly 
on T, thus making the oscillation a very reliable temper- 
ature probe in the low temperature regime where current 
techniques fail. 



VII. SUMMARY AND OUTLOOK 

In this article, we have investigated the phenomenon of 
quantum oscillations for a non-interacting Fermi gas sub- 
ject to rotation or to an artificial magnetic field. We have 
considered the angular momentum of the gas, which is di- 
rectly analogous to the magnetization in the solid-state 
context, and hence for which de Haas - van Alphen oscil- 
lations are expected as a function of e.g. atom number 
or chemical potential. 

We have uncovered important differences between 
atomic gases and the solid-state context, especially in 
the limit of fast rotation a; ^ cjq for a two-dimensional 
system. These differences are due to the different de- 
generacies of the effective Landau levels. As a result, 
the pronounced periodic jumps in the magnetization per 
electron of a two-dimensional electron gas are replaced in 



the cold fermionic gas context by milder cusp-like singu- 
larities of the angular momentum versus atom number. 

In contrast, we have shown that pronounced oscilla- 
tions do show up in the opposite limit of slow rotation. 
The angular momentum is an observable which is par- 
ticularly prone to oscillations, in contrast to the mean 
square radius of the gas for example (related to the rigid- 
body moment of inertia) , which only displays small shell- 
filling effects for very small systems. We performed a de- 
tailed investigation of the dependence of these oscillations 
on experimentally relevant parameters. We concluded 
that they are observable up to rather large atom num- 
bers (several thousands) provided the harmonic confine- 
ment along the rotation axis is weak. The temperature 
must also be low compared to hw^/ks — 1t-{ijJq — a;)/fcs, 
so that a significant effect for reasonably large systems 
is typically found only for T/Tp < 0.1. In the solid- 
state context, quantum oscillations are a tool of choice 
for determining the Fermi surface of clean-enough mate- 
rials [16, 33[. Such use is also possible in the context of 
atomic gases, but the technique can probably not com- 
pete with adiabatic band mapping [34[ . Quantum oscilla- 
tions could, however, provide a useful method for in-situ 
thermometry in the regime of low T /Tp where the ampli- 
tude of the oscillations depends strongly on temperature 
and where time-of-fiight or in-situ measurements of the 
density profile are less sensitive. 

We proposed two experimental protocols which could 
be used to observe these quantum oscillations experimen- 
tally. These protocols elaborate on previous propositions 
to measure the angular momentum of the gas [7, 30[ in- 
spired from the historic Andronikashviili experiment [35[ 
on superfiuid helium. There are several extensions of 
the present work that are worth considering. One is the 
effect of interactions on the quantum oscillations. In- 
teractions and the associated finite lifetime of the quasi- 
particles will lead to a damping of the oscillations which 
will reduce their visibility [23[. At a phenomenological 
level, this can be accounted for by a Dingle-like factor, 
which in our case reads e~''''"+/"++"-/"-^ in the oscillat- 
ing contributions (assuming that the broadening due to 
interactions is the same for all energy levels). The char- 
acteristic damping 7 can be interpreted as the inverse 
lifetime of quasiparticles and its measurement (e.g. as a 
function of scattering length) could thus provide useful 
information on this lifetime. 

Another interesting direction is to investigate quantum 
oscillations in the presence of an optical lattice and an 
artificial gauge field, exploring in particular the possible 
signature of non-trivial topological states. 
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Appendix A: Poisson sums 

Poisson summation is useful to compute sums of the 
form 



oo 

n— no 



(Al) 



and has been used by Landau to predict quantum oscil- 
lations for clean electronic systems [1]. To describe the 
results presented in the main text, a sum over three vari- 
ables n^,n_,nz is needed. In full generality, the result 
of this procedure is given by : 



[0,+ooP 



2^ E 



du(j){u)e 



2i7rk-u 



, (A2) 



where h' = {n+,n-,nz) in the left hand side of eq. A2 
and k — (k^,k-,kz) with k±,kz positive integers. The 
expression (A2) already displays the separation between 
monotonous and oscillatory contributions. 

Let us apply this procedure to the average of a generic 
quantity $ which reads 



m = 



(A3) 



Using equation (A2), this average can be separated 
into two contributions 



Here the monotonous part is given by 

$(u) 



du ■ 



'[0,+oop 1 +exp [/3(ea - pi)] 
and the oscillating contribution by 



(A4) 



(A5) 



.MO 



k, fc _ J.- \ 



(A6) 



In the last expression, the quantity T{k) is given by : 



+00 f+oo 

de I dy' 



-t-oo 



2iir 



1 -I- e/3(£+eo) 

e + ii + y'/'S + z'/9 e + n - y' /'S + z' /9 e + ^j.-2z'/9 



hWz/i 



, (A7) 



r 



where the indices k± , kz are all positive integers and 
(fc+,fc_,fc,)^ (0,0,0). 

This analysis shows how the oscillating behavior arises in 
thermodynamical quantities such as the angular momen- 
tum studied here. Furthermore, one sees that the tem- 
perature dependence of the oscillating part is contained 
in J- , which is the convolution between the Fourier trans- 
form of $ and the one of the occupation number, which 



behaves like ( sinh 



keT 



Thus, one can expect an exponential damping of the 
oscillating term due to temperature. 



lating part are quite complicated, and not very illumi- 
nating. Nevertheless, the decomposition (A4) provides 
rather simple expressions for the non oscillating contri- 
butions. 



This computation also shows that the period of the os- 
cillations are related to the inverse frequencies and 



jri-. In our case, the axial trapping frequency is usu- 
ally quite small, making the corresponding period very 
long, and thus not observable in realistic temperature 
conditions. In most cases, the expressions for the oscil- 



Appendix B: Exact expressions for non oscillating 
contributions 



The generalization of (A2) for several variables leads 
to the following expressions for the non oscillating con- 
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tributions. at any temperature : 
1 



N 



-J 



1 



21 



21 + I3huja{l - ol^)J 

- - I du ^ = 

6 io 1 + e"-« 



/3(A'-£o)1 



-Li4[— 6' 



-Li3[-e'5(^-=«)] 



with Li„ the n-th polylogarithm function[19]. 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 



These expressions show in particular that at large /i, 
the populations of the radial harmonic oscillators are 
such that = TT^- Also, we see that the devia- 

tion of the ratio between the angular momentum of the 
gas and its classical expression goes like ^, thus like 
jY-1/3 j3gj macroscopic samples. The strictly zero 
temperature result for the number of particles leads to 
Ef = /k^o(6a^7V)i/3 [36], as expected. 



Appendix C: Angular momentum of a classical gas 

We consider a gas of classical particles interacting 
through the potential V trapped in a potential Vtrap- 
The whole gas is rotating at constant angular frequency 
Lo around the z-axis : 

^ ^ E + ^^rap in )-L,,^oj + J2v{\n-r,\). (CI) 

i ij 

The average angular momentum of the gas is given by : 

5Tr (e-'3«) 



duj 



(C2) 



One can rewrite the Hamiltonian (CI) as : 



^ = 2^ ^ + Vtrapin) - -Muj {x +y) 



(C3) 



showing that, analogously to Bohr- van Leeuwen theorem 
(see eg [37]), with an appropriate change of variables in 
the partition function, one can immediately prove that 
the angular momentum is given by eq. (10). Further- 
more, one also sees that the interaction potential V does 
not play any role in deriving the expression of the angu- 
lar momentum, showing that the result (10) also holds 
for interacting particles. 



[1] L. Landau, Z. Phys 64, 629 (1930). 

[2] W. de Haas and P. van Alphen, Comm. Phys. Lab. Leiden 
212a (1930). 

[3] N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 

J.-B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 

and L. Taillefer, Nature 447, 565 (2007). 
[4] A. L. Fetter, Rev. Mod. Phys. 81, 647 (2009). 
[5] N. Cooper, Advances in Physics 57, 539 (2008). 
[6] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[7] F. Chevy, K. W. Madison, and J. Dahbard, Phys. Rev. 

Lett. 85, 2223 (2000). 
[8] M. Zwierlein, J. Abo-Shaeer, A. Schirotzek, C. Schunck, 

and W. Ketterle, Nature 435, 1047 (2005). 
[9] Y.-J. Lin, R. L. Compton, K. Jimenez-Garcia, J. V. 

Porto, and I. B. Spielman, Nature 462, 628 (2010). 
[10] J. Dalibard, F. Gerbier, G. JuzeUunas, and P. Ohberg, 

Rev. Mod. Phys. 83, 1523 (2011). 
[11] S. Stringari, Phys. Rev. Lett. 76, 1405 (1996). 
[12] F. Zambelh and S. Stringari, Phys. Rev. Lett. 81, 1754 

(1998). 

[13] P. C. Haljan, B. P. Anderson, L Coddington, and E. A. 

Cornell, Phys. Rev. Lett. 86, 2922 (2001). 
]14] A. J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 
]15] A. J. Leggett, Journal of Statistical Physics 93, 927 

(1998). 



[16] A. Abrikosov, Introduction to the theory of normal met- 
als, SoHd State Physics: Supplement (Academic Press, 
1972). 

[17] C. Cohen-Tannoudji, B. Diu, F. Laloe, and B. Diu, 
Quantum Mechanics (2 vol. set) (Wiley-Interscience, 
2006). 

[18] This monotonous deviation is also thought to be rather 
a quantum effect than a finite size one, since the classical 
expression for the angular momentum remains true even 
for a single classical particle. 

[19] M. Abramovitz and L Stegun, Handbook of mathematical 
functions, 10th edition, (Dover, 1972). 

[20] J. Schneider and H. Wallis, Phys. Rev. A 57, 1253 (1998). 

[21] T.-L. Ho and C. V. Ciobanu, Phys. Rev. Lett. 85, 4648 
(2000). 

[22| P. Vignolo and A. Minguzzi, Phys. Rev. A 67, 053601 
(2003). 

[23| G. M. Bruun and K. Burnett, Phys. Rev. A 58, 2427 
(1998). 

[24[ D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 
(1997). 

[25| C. Weisbuch and B. Vinter, Quantum Semiconductor 
Structures: Fundamentals and Applications (Academic 
Press, 1991). 

[26[ S. Nascimbene, N. Navon, K. Jiang, F. Chevy, and C. Sa- 
lomon, Nature 463, 1057 (2010). 



14 



[27] N. Navon, S. Nascimbene, F. Chevy, and C. Salomon, 
Science 328, 729 (2010). 

[28] V. Bretin, S. Stock, Y. Seurin, and J. Dalibard, Phys. 
Rev. Lett. 92, 050403 (2004). 

129] K. W. Madison, F. Chevy, W. Wohlleben, and J. Dal- 
ibard, Phys. Rev. Lett. 84, 806 (2000). 

130] N. R. Cooper and Z. Hadzibabic, Phys. Rev. Lett. 104, 
030401 (2010). 

131] S. T. John, Z. Hadzibabic, and N. R. Cooper, Phys. Rev. 
A 83, 023610 (2011). 



[32] J. Dalibard, (2000), unpublished. 

[33] N. Ashcroft and N. Mermin, Solid State Physics (Saun- 
ders College, Philadelphia, 1976). 

[34] M. Kohl, H. Moritz, T. Stoferle, K. Giinter, and 
T. Esslinger, Phys. Rev. Lett. 94, 080403 (2005). 

[35| E. AndronikashviU, J. Phys. USSR 10, 201 (1946). 

[36[ C. Pethick and H. Smith, Bose-Emstein condensation in 
dilute gases, (Cambridge University Press, 2008). 

[37| S. Blundell, Magnetism in condensed matter, Oxford 
master (Oxford University Press, 2001). 



